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ABSTRACT 

Theories of high energy radiation production in quasar jets can be verified by 
studies of both time- averaged spectra and variabihty patterns. While the former 
has been explored extensively, the latter is in its infancy. In this paper, we 
study the production of short-term flares in the shock-in-jet model. We examine 
how the flares' proflles depend on such parameters as shock/ dissipation lifetime, 
electron-injection time profile, adiabaticity, and half-opening angle of the jet. In 
particular, we demonstrate the large difference between fiare profiles produced in 
the radiative and adiabatic regimes. We apply our model to the ~day timescale 
fiares observed in optically violently variable (OVV) quasars, checking whether 
the external-radiation-Compton (ERC) model for 7-ray fiares at energies > 30 
MeV (EGRET range) can be reconciled with the fiares observed at lower energies. 
Specifically, we show that the strict correlation between X-ray and 7-ray fiares 
strongly supports the dominance of the synchrotron self-Compton mechanism 
in the X-ray band. We also derive conditions that must be satisfied by the 
ERC model in order to explain a lag of the 7-ray peak behind the optical one, 
as claimed to be observed in PKS 1406-076. Finally, we predict that in ERC 
models where the MeV peak is related to the break in electron distribution due 
to inefficient cooling of electrons below a certain energy, the fiares should decay 
significantly more slowly in the soft 7-ray band than at energies greater than 30 
MeV. 

Subject headings: galaxies: quasars: general — galaxies: jets — radiation 
mechanisms: nonthermal — gamma rays: theory — X-rays: general 



-3- 

1. INTRODUCTION 

The blazar phenomenon is well estabhshed to be related to nonthermal processes 
taking place in relativistic jets on parsec/subparsec scales and viewed by observers located 
within or nearby the Doppler cone of the beamed radiation. Thus, studies of the properties 
of blazar radiation, such as spectrum shape and variability, provide exceptional tools 
for exploring the deepest parts of extragalactic jets, their structure, physics and origin. 
However, in order to accomplish this, one needs first to identify the dominant process 
responsible for the production of 7-rays. 

The 7-rays, with a luminosity peak located at photon energies > 1 MeV, form the 
high energy spectral component and, together with the synchrotron component peaking in 
the IR-X-ray range, constitute the characteristic two-component spectrum of blazars (von 
Montigny et al. 1995; Fossati et al. 1998). The 7-ray fluxes show rapid and high-amphtude 
variability, and in many OVV and HP (highly polarized) quasars, at least during their high 
states, they reach luminosities 10-100 times larger than in the lower energy spectral bands. 

There are several radiation mechanisms that can contribute to the production of 7-rays 
in relativistic jets. One, suggested by Konigl (1981), is the synchrotron self-Compton (SSC) 
process. Others involve Comptonization of external radiation fields. The candidate fields 
are: direct radiation from the disk (Dcrmcr & Schlickeiser 1993); broad emission lines 
(BELs) and near-IR radiation of hot dust (Sikora, Begelman & Rees 1994; Blazejowski et 
al. 2000 [hereafter, B2000]); rescattered central X-ray radiation (Blandford & Levinson 
1995); and externally rescattered/reprocessed synchrotron radiation of the jet (Ghisellini & 
Madau 1996). Gamma-rays can also be produced by synchrotron radiation of an extremely 
relativistic population of electrons/positrons, with the maximum random Lorentz factor 
required to be at least 10^ — 10^. Such a population of electrons/positrons can be the 
product of a synchrotron pair cascade, induced by > 10^ GeV protons via the photo-meson 
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process (Mannheim and Biermann 1992). In some circumstances, a certain contribution 
to the 7-ray band can also be provided by synchrotron radiation of muons (Rachen & 
Meszaros 1998; Rachen 2000) and protons (Aharonian 2000; Miicke & Protheroe 2000). 

Hereafter we focus on quasars only. In these objects the energy density of the external 
diffuse radiation field is very high and, as amplified by in the jet comoving frame, where 
r is the bulk Lorentz factor of radiating plasma, can easily dominate over the energy 
density of magnetic fields carried by the jet and the energy density of synchrotron radiation 
produced in the jet (see, e.g., Sikora 1997). Thus, the radiative output of quasars is very 
likely to be dominated by 7-rays produced by Comptonization of external diffuse radiation 
fields. This so-called external radiation Compton (ERC) model naturally explains the 
location of the 7-ray peak/break, which in OW and HP quasars is observed to be in the 
1-30 MeV range. The MeV break, according to the one-zone ERC model version proposed 
by Sikora et al. (1994), simply matches the break in the energy distribution of electrons 
caused by their inefficient cooling at lower energies. The location of this break requires the 
flares to be produced ~ 0.1 — 1.0 pc from the core, a distance which corresponds nicely with 
the observed flare time scales ~ days, provided that the formula r ~ ^^ctfi applies, where r 
is the distance from the central engine and tfi is the time scale of the flare. 

The specific distance range of 7-ray production can be related to the distance range 
over which the collision of two inhomogeneities moving with different velocities is completed. 
According to this scenario, the variability of blazars is modulated by the central engine via 
instabilities in the innermost parts of an accretion disk or magnetic eruptions in the corona 
(Sikora & Madejski 2000). If separations and lengths of inhomogeneities are of the same 
order, then the number of reverse-forward shock pairs enclosed within a dissipative zone is 
about and the number of shocks observed at a given moment is of order 1 (Sikora et al. 
1997). This can explain high-amplitude fluctuations of blazar light curves, with occasional 
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flares exceeding the background blazar radiation by more than a factor 3. 

In this paper, we use the shock-in-jet model to check whether 7-ray flares produced 
by the ERC process can be reconciled with those observed at lower energies (§3). In 
§3.1, we show that the similar rates at which 7-ray and X-ray flares decay, as observed 
in 3C279 (Wehrle et al. 1998), require the production of X-rays to be dominated by the 
SSC process, as was suggested by Inoue and Takahara (1996) and by Kubo et al. (1998). 
In §3.2, we study the possible time-lags of flares produced at different frequencies and 
discuss conditions required to reproduce the optical-7-ray lag claimed to be recorded in 
PKS 1406-076 (Wagner et al. 1995). Finally, in §3.3 we make predictions regarding the 
variability of 7-rays at energies > 30 MeV vs. variability in the soft 7-ray band, around 1 
MeV. These predictions are aimed at verifying whether our interpretation of the MeV break 
in terms of electron cooling effects is correct. These studies are preceded by analyses of 
how the flare profile depends on such model parameters as shock lifetime, particle injection 
function, adiabaticity and jet opening angle (§2). Our results are summarized in §4. 

2. FEATURES OF THE MODEL 

2.1. Model cissumptions 

To study flares in blazars we adopt the shock-in-jet scenario, in which individual flares 
are produced by shocks formed due to velocity irregularities in the beam and traveling 
down the jet with relativistic speeds. Our model assumptions/ approximations are (here, 
as elsewhere in the paper, primed quantities denote measurements made in the co-moving 

frame of the source): 

— nonthermal plasma producing flares is enclosed within thin shells, having a radial 
comoving width. A', much smaller than their cross-sectional radius a; 
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— shells propagate down the conical jet with a constant Lorentz factor F; 

— magnetic fields, carried by the beam, scale with distance as oc 1/r; 

— both magnetic field intensity and particle distribution are uniform across the shell; 

— relativistic electrons/positrons are injected into the shell within a finite distance range, 

^finj-i which is equal to c/Stcoiu where /Stcoii is a timescale of the collision; 

— the injection rate is parameterized by Q = Kj'^, for 7m < 7 < Imax, and Q oc 7"^, for 

7 < 7m) where 7 is the random Lorentz factor of an electron/positron; 

— radiative energy losses of relativistic electrons/positrons are dominated by Comptonization 

of external radiation fields, synchrotron radiation, and Comptonized synchrotron 
radiation; 

— the observer is located at an angle 9obs = 1/r from the jet axis. 

The above assumptions can be justified as follows. A shell in our model approximates 
the geometrical site of relativistic electrons and positrons. Such a site can be identified 
with the space between forward and reverse shocks, formed due to the collision of two 
inhomogeneities moving down the jet with different velocities. The speed of the shell should 
then be identified with the speed of the contact surface, while the width of the shell is just 
the distance between forward and reverse shock fronts. This distance is of course a function 
of time and in the contact surface frame is given by the formula A' = + (3'j.)cAt', where 
/S'fC and /3^c are forward and reverse shock velocities, respectively, and At' is the time since 
the collision started. Comparing this with the cross-sectional radius of the shell, a = r/T, 
one finds that 

-^(P'f+P'r)^- (1) 

a r 

As was demonstrated by Komissarov and Falle (1997), for intrinsically identical 
inhomogeneities (the same density and pressure) having Lorentz factors differing by less 
than a factor two, the shock velocities as measured in the contact frame are non-relativistic 
and, therefore, our approximation A' -C a is justified, provided Ar^^j < r. 
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Another geometrical assumption concerns the cone geometry. As VLBI observations 
show, parsec scale radio jets have a tendency to be wider at subparsec distances than on 
larger scales (Lobanov 1998) and, within a finite distance range, they are much better 
approximated by conical geometry than by cylindrical geometry. We should note that even 
for H^Tinj/r <^ 1 it is important to take into account the radial divergence of the jet, simply 
because radiation from shell elements moving in different directions with respect to the 
observer are differently Doppler boosted. (Of course, given possible departures of the front 
and rear surfaces of the colliding inhomogeneities from a spherical shape, one can expect 
that in general the shock surfaces will not be exactly spherical. However, if these departures 
are 5A' <C Arj„j/r, the effect on fiare profiles is negligible and the sphere approximation 
can be justified.) 

In all of our models except those presented in Figure 2d the half-opening angle of 
the jet is taken to be 9j = l/F. This choice is somehow arbitrary, but our results are not 
expected to be dramatically different if the real scaling is larger or smaller by a factor 
three. There are several more or less direct arguments that 9j ~ l/F. The deepest radio 
observations of a quasar jet are those of Cygnus A, a radio galaxy that hides in its center a 
powerful quasar (Antonucci, Hurt, & Kinney 1994). The jet opening angle in this object on 
0.1 parsec scale is about 7 degrees (Lobanov 1998; Krichbaum et al. 1998). This is equal 
to the full Doppler angle, 2/F, if F ~ 17, whereas bulk Lorentz factors of 7-ray quasars, as 
deduced from ERG models, are enclosed within the range 10 — 20 (sec, e.g., Ghisellini et 
al. 1998). For 9j -C l/F, the SSC radiation, calculated self-consistently assuming an ERG 
origin of 7-rays, would exceed the observed soft X-rays; while 9j ^ l/F would increase jet 
energetic requirements up to a level difficult to reconcile with theories of the central engine. 
Furthermore, with 9j ^ l/F the number of blazars would be too large (Padovani & Urry 
1992; Maraschi & Rovetti 1994). 

We should also comment on our assumptions regarding the particle acceleration 
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process. One is that the injection of relativistic electrons is described by a two-power-law 
function, with a break at 7^- This break corresponds roughly with the average energy of 
injected electrons, and the low-energy tail below the break mimics a limited efficiency of 
the electron pre-heating process. It should be noted that the model output parameters are 
insensitive to the exact slope of this very hard low-energy tail, and that observationally 
it can be imprinted only in the hard X-ray band. In the synchrotron component it is 
invisible because of the synchrotron-self- absorption process, and in soft-to-mid-energy 
X-rays it is obscured by the much stronger SSC component. Another assumption, which 
is not explicitly listed above, concerns the acceleration time scale. Since the time scale of 
particle acceleration cannot be longer than the time scale of radiative losses of even the 
most relativistic electrons, and since the latter is known from modeling electromagnetic 
spectra of OW/HP quasars to be ^ a/c, the particle acceleration can be approximated as 
an instantaneous process, represented in the electron kinetic equation as a separate term 
called the injection function and denoted by Q (B2000). The particle acceleration time 
scale should not be confused with the time scale of particle injection, a term which we 
reserve hereafter for the time over which the shock operates, and during which particles are 
accelerated continuously. 

2.2. Evolution of the electron energy distribution 

Evolution of electrons is described by a continuity equation (Moderski, Sikora & Buhk 
2000; B2000) 




(2) 



where the rate of electron/positron energy losses is 




(3) 
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P — -\/r^^^/r, and dr — PcTdt'. The second term on the right-hand side of eq. (3) 
represents the adiabatic energy losses, with A — 1 ior 3D expansion and A = 2/3 for 2D 
expansion. 

In Figure 1 we show examples of electron evolution for the case of radiation energy 
losses {dj/dt')rad 7^. Plotted is N-yj"^ vs. 7/70, where A^^ is the number of electrons per 
unit of energy and 7,. is the energy of electrons for which the radiative cooling time scale is 
equal to the injection time scale. Noting that the electron radiative energy loss time scales 
are given by 

trad 



7 



{dl/dt') 



r (4) 

rad 



where \d^/dt'\j.^^ — Irf'^, one can find that 



For b dependent on r, 7c is computed for b taken at r = Tq + Arj„j. 

As one can see in Figure 1, for r < Arj„j the number of electrons is increasing, being 
saturated first at high energies, then at lower energies. The slope of the saturated radiative 
part of the electron energy distribution is given by ~ 7~^~^, while at energies 7 < 7c 
the electron distribution follows the injection function. After injection stops, the energy 
distribution of electrons above the break steepens and the break itself moves down to ~ 7c, 
or even lower energies, depending on whether adiabatic losses are taken into account or not 
(compare the left panels with the right panels). 

We also illustrate in Figure 1 how the electron evolution is different depending on 
whether the electron radiative energy losses are constant or drop with distance. In the 
latter case, the high energy parts of the electron distribution saturate initially at lower 
amplitudes, due to the fact that the energy losses are faster at the beginning of the injection 
process than they are later. This amplitude increases with time, reaching a maximum at 
the end of the injection process. 
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Finally, we note that the two cases considered here, one with d'y/dt' — const and one 
with d^jdt' oc 1/r^, can be related to two possible configurations for the source of radiation 
surrounding the jet. In the first case the radiation energy density in the comoving frame of 
the shell is dominated by external radiation sources (BEL-clouds or hot dust) located at a 
larger distance from the central engine than the shell, while in the second case both are 
at a similar distance and the energy density is scaled by ^L/r^, with the fraction ^ of the 
reprocessed/rescattered central luminosity, L, being constant. 

2.3. Flares 

Using the algorithm developed by B2000 to compute the observed radiation from 
relativistically and "conically" propagating shells, we analyze in this section the dependence 
of fiare profiles on: electron injection distance range, Arj„j/ro; injection time profile, 
adiabaticity, A\ and jet opening angle, Qy Flares, presented in Figure 2, are produced by 
the ERC process and are computed under the assumption that electron radiative losses 
are strongly dominated by just this process. All panels except (b) show flares which 
are produced in the radiative regime, i.e., at an observed frequency that is contributed 
by electrons with 7 ^ 7c. One can easily check that "radiative" flares reach maximum 
at a time when the observer is receiving the flrst photon from a radius ro + Arj„j. For 
^Tinj/vQ = 1, this is at t/to = 1, where to = ^o/(l + P)PcT'^ — ro/2cr^. The small shifts of 
the peaks to the right of this value are caused by the flnite time of electron coohng. 

As one can see from panel (a), flares produced by long lasting injections, i.e., with 
Arj„j > ro, have very shallow peaks. Such flares are never observed in blazars, suggesting 
that Arj„j < Tq. On the other hand, those with very small l\rinj/rQ are predicted to be 
more asymmetric than observed flares, showing much faster increase than decay, the latter 
determined by the light travel time effect related to the transverse size of the shell. 

Panel (b) shows "adiabatic" flares, produced by electrons with 7 <^ 7^. One can see 
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that even those with maximal adiabatic losses (case A — 1) are decaying significantly slower 
than radiative flares. This feature can be used to verify X-ray production mechanisms (see 
§3.1). 

As panel (c) shows, the dependence of flare proflles on the time proflle of electron 
injection, Q{t'), is rather weak. This is because this dependence is diluted by the 
transverse- size light travel time effect. Only for Arj„j/ro > 1, for very narrow jets (see 
panel [d]) or for the observer located far from the jet axis can one try to recover the 
injection history Q{t') from the observed fight curves. To better illustrate this effect we 
present in Figure 3 fiares produced by a double-triangle injection profile. As one can see, 
the double-peak structure of the flare becomes more diluted as the injection time becomes 
smaller or the observer moves out from the jet axis. 

2.4. Frequency-dependent lags 

Very useful constraints on jet structure and radiation models can be provided in the 
future by detailed observations of time lags between flares detected in different spectral 
bands. In this section, we demonstrate the predictions of our model regarding the lags: at 
different frequencies within the synchrotron component; at different frequencies within the 
ERG component; and between synchrotron, ERG and SSG components. 

Our model parameters are chosen to match flducial spectra of 7-ray quasars during 
their high states (see, e.g., Fossati et al. 1998). The spectrum evolution and lags predicted 
by our model are presented in Figure 4. The observed frequencies of flare production in the 
synchrotron component (panel a) and in the ERG component (panel b) were chosen to be 
contributed by the same electron energies at the beginning of the flare. At later moments 
this is no longer true. This is because the synchrotron critical frequency is proportional 
to the magnetic fleld intensity, Ugyn oc T'j'^B', thus the synchrotron spectrum produced 
by electrons from a given range of the energy distribution moves to the left. In the ERG 
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case, where uerc ~ ^'^^'^i^ext and Uext ~ const (~ 10 eV for BELs and ~ 0.3 eV for near 
IR radiation of hot dust), the spectrum doesn't drift with time. This is the main reason 
why synchrotron flares observed at different frequencies are predicted to show lags (panel 
a), while the effect is negligible in the ERC component (panel b): because the critical 
synchrotron frequency is dropping, it takes longer to build up the population of electrons 
necessary to produce a peak at higher frequencies. Of course, one can see that the ERC 
flare profiles also depend on frequency, but this is caused by the energy dependence of the 
electron cooling rate, and the effect is significant only at frequencies at which radiation is 
produced by electrons with 7 < 7c. 

There are two additional factors that differentiate the profiles of synchrotron 
flares from those of ERC flares. The first is that the synchrotron cooling rate is 
{d^/dt')syn oc B'^ oc 1/r^, whereas (d'y/dt')ERC is assumed in this model to be constant. 
Since r is increasing with time as the shock propagates, this makes the synchrotron flare 
peak earlier. The second is that the synchrotron radiation in the "comoving" frame is 
isotropic, whereas the ERC radiation is anisotropic (see Dermer 1995; Sikora 1997; B2000). 
Direct comparison of flares produced by different processes is provided in panel (d), at 
frequencies which are marked by arrows in panel (c). 

3. CONFRONTATION WITH OBSERVATIONS 

3.1. X-rays vs. 7-rays 

Unfortunately there were only a few multi-wavelength campaigns, during the lifetime 
of the EGRET/CGRO detector, that permitted measurements of time lags or correlations 
between flares observed in the 7-ray band and at longer wavelengths. Undoubtedly, the 
most successful was the one at the beginning of 1996, aimed at 3C 279 (Wehrle et al. 1998). 
During this campaign, simultaneous coverage was obtained in practically all frequency 
bands. 
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A particularly interesting result of the campaign was the very close match between 
the X-ray and 7-ray light curves during the Feb-96 outburst event. There was no 
time lag between the peaks, and both decayed at the same rate (Lawson, McHardy & 
Marscher 1999). This is not what one would expect, assuming that ERC dominates the 
production of radiation down to the X-ray bands, because the production of X-rays by 
ERC involves low energy electrons, which radiate in the adiabatic regime. As was shown in 
Figure 2b, adiabatic flares decay significantly more slowly than radiative ones. Therefore, 
the observed 7-X-ray correlation seems to exclude the production of X-rays by the ERC 
process, unless the decay rate is determined by the deceleration of the shell and/or a change 
in its direction of motion. 

We applied our ERC model for 7-ray production during the Feb-96 outburst and 
found that the main contribution to the X-ray band is coming from the SSC process. This 
interpretation of the origin of the X-rays is consistent with both the time-averaged fit to 
the spectrum and the correlation between the 7-ray and X-ray flares (see Figure 5 upper 
panels). It should be noted, however, that our model predicts comparable X-ray and 7-ray 
flare amplitudes, while observations show that the 7-ray flux jumped during the outburst 
by a larger factor than the X-ray flux did. But this can be explained by dilution of the 
X-ray flare by radiation produced at larger distances in the jet. It is worth mentioning that 
even stronger dilution is required in the optical range, where the outburst was hardly visible 
(Wehrle et al. 1998). 

In Figure 5c we also present our X-ray flare proflle superimposed on observations taken 
during the Feb-96 campaign. Data points represent 2 — 10 keV flux measured by RXTE 
(Lawson et al. 1999). As one can see the "flt" is not perfect, especially on its rising side, 
and shows signiflcant deviation from the exponential symmetrical proflle used by Lawson 
et al. (1999). We attribute this discrepancy to two factors. First of all, the model flare in 
Figure 5c is produced by a rectangular injection proflle. As can be seen from Figure 2c, 
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triangular injection profiles tend to produce more symmetrical flares. Secondly, the rise of 
the X-ray flux before the 7-ray flare may indicate the existence of a flare precursor. Such 
precursors may be produced by Comptonization of external radiation by cold electrons 
necessarily present in inhomogenities before their collisions and shock formation (Sikora et 
al., in preparation). 

3.2. Optical— 7-ray lag: A mirror model? 

Another noticeable multiwavelength event was recorded in PKS 1406-076 (Wagner et 
al. 1995). In this case, unlike in the Feb-96 flare in 3C279, the 7-ray flare was clearly 
accompanied by an optical one. Detailed analyses of this event suggest that the peak of 
the 7-ray flare lagged the optical peak by about one day. If this is true and this type of 
lag is confirmed by future observations, it would have dramatic implications for radiation 
scenarios in blazars. Whereas our model predicts an optical-7-ray lag (see panel [d] in 
Figure 4), it cannot reproduce such a pronounced lag as claimed for PKS 1406-076 during 
the Jan- 93 outburst. 

A model which can be reconciled with the observed lag is the "mirror" model, suggested 
by GhiseUini & Madau (1996). According to this model the synchrotron and Compton 
flares are produced by the same source (i.e., the same shock), but at different distances. 
Synchrotron flares are produced closer to the center. Radiation from these flares is 
scattered/reprocessed by BEL clouds, then Comptonized by the propagating shock, giving 
rise to the delayed 7-ray flares. In modeling such a scenario, it is necessary to note that the 
energetics of nonthermal flares produced by electrons in the radiative regime is determined 
by the rate of relativistic electron injection in the jet, rather than by the intensity of the 
external radiation held. Also, the model should take into account the fact that the 7-ray 
flux dominated the radiation output of PKS 1406-076 during the entire Jan-93 event. 
Hence, the 7-ray flare must be modulated by Q{t'), whereas the synchrotron flare proflle 
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depends on both Q{t') and the time dependence of u'^/u'^^f^, where u'^ is the energy density 
of the magnetic field in the propagating source, and w^i// ^^e energy density of externally 
rescattered/reprocessed synchrotron radiation of the jet. Our model results are presented 
in Figure 6. 

We would also like to comment on the possible nature of the material that scatters and 
reprocesses synchrotron radiation outside the jet. Can it really consist of clouds producing 
broad emission lines? Since the flux of synchrotron radiation near the jet is much larger 
than the flux of radiation from the central engine at this distance, the clouds near the 
jet are expected to be ionized to a much higher level than those which are not exposed 
to jet radiation. The question is then whether they can produce strong emission lines. 
Furthermore, there are some arguments in favor of the BEL region being identified with 
accretion disk winds (Emmering, Blandford & Shlosman 1992; Murray & Chiang 1997; 
Proga, Stone & Kallman 2000; Nicastro 2000). If this is the case, then explaining the 
optical-7-ray lag through a mirror model would require a different reprocessor/rescaterrer 
than BEL clouds. Alternatively, the increase of u'^/u'^^^^ required between the regions 
producing the synchrotron and the Compton peaks can be obtained by posting a source of 
seed photons in the vicinity of the jet, at a distance beyond the region where the shock is 
launched. Such a source could be provided by a supernova explosion, but the probability of 
supernovae on parsec scales in the vicinity of the jet seems to be rather low. 



3.3. On the origin of the MeV break 

In all previous sections it was assumed that the break in the electron injection function 
(7^) is located at lower energies than the break caused by inefficient radiative cooling of 
electrons with 7 < 7c. With this assumption, the characteristic break of the high energy 
spectra in the 1-30 MeV range is related to the break of the electron energy distribution at 
7c. As was discussed by Sikora (1997), such an interpretation of the MeV break is consistent 
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with the observed flare time scales. Assuming that the distance of flare production is of 
the order Vfi ^ ctfiV^, one can find that flares lasting ~days are produced at distances 
10^^ — lO^^cm, and that at such distances ERC radiation is inefficient ai hu < 1 — 30 MeV 
(see also Sikora et al. 1996 for discussion on coohng efficiency in different regions of the jet). 

However, the time resolution of recent high-energy experiments is not good enough 
to reject the possibility that the observed flares are superpositions of several short-lasting 
flares, which are produced much closer to the central engine than the ~ 1-day variability 
time scales suggest. Due to the stronger magnetic and external radiation fields at smaller 
distances, the value of 7c would be lower and the break would be imprinted at hv <^ I 
MeV. In order to explain the observed MeV breaks, one would then have to assume that 
this break is related to the break in the electron injection function. 

There are two ways to distinguish between these possibilities, by studying spectral 
slopes and by studying flare profiles, both below and above the MeV break. The first 
approach is based on the fact that for 7^ < 7c, the slope of the radiation spectrum produced 
by electrons with 7m <^ 7 <S 7c should be harder by Sa — 0.5 than the slope of the radiation 
spectrum produced by electrons with 7 3> 7c. But in the case 7c < 7m, the slope of the 
radiation spectrum produced by electrons with 7c <C 7 -C 7m is predicted to be a = 0.5, 
independent of the spectral slope produced at u ^ ^(7m) (see, e.g.. Sari, Piran & Narayan 
1998). The second approach, proposed in this paper, uses the fact that adiabatic flares 
decay much more slowly than radiative ones (sec Figure 2). Thus, in the case 7^ < 7c, the 
flares observed in OVV/HP quasars at < MeV energies should decay signiflcantly more 
slowly than in the case 7c < 7m- We present our predictions for two models, one in which 
d'f/dt' = const (Figure 7) and one in which d^/dt' oc (Figure 8). In the latter model 
the difference between the two cases is narrowed by the fact that 7c is not constant. 
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4. CONCLUSIONS 

We have studied flare production by thin shells propagating at relativistic speeds down 
a jet and applied this model to 7-ray quasars. Our main results can be summarized as 
follows: 

- the sharp nature of flares observed in OVV/HP quasars (sec, e.g., Mattox et al. 1997; 

Wagner et al. 1995; Wehrle et al. 1998) seems to conflict with models in which the flare 
production time scale is much longer than the light travel time across the source. Such 
models produce flares which have profiles determined by the injection history Q{t'), 
and, unless Q is a very sharply peaked function of time, these flares are too shallow to 
be consistent with observations; 

- the fast decay of the X-ray flare observed during the Feb-96 event in 3C 279 suggests that 

the production of X-rays in this object is dominated by the SSC process; 

- the claimed lag of the 7-ray flare behind the optical flare during the Jan-93 outburst 

in PKS 1406-076 requires a scenario in which both the energy density of the external 
radiation field and the injection rate of electrons increase between the production of the 
synchrotron peak and of the ERC peak; 

- interpretation of the MeV break in terms of inefficient cooling of electrons radiating at 

lower energies predicts a much slower decay of MeV flares than of GeV flares. 
The above conclusions can be weakened, or even invalidated, if the flare decay profiles are 
determined by deceleration of the shell or a changing direction of motion. Such a possibility 
can be tested by future simultaneous observations of flares at GeV, MeV and keV energies. 
Since kinematically determined flare decays should be achromatic, finding an example in 
which X-ray flares decay as quickly as GeV flares, whereas MeV flares do not, would argue 
against the role of kinematics in shaping the flare profiles. 



This project was partially supported by Polish KBN grant 2P03D 00415, NSF 



-18- 



grants AST-9876887 and AST-9529175, and NASA grant NAG5-6337. M.S. thanks JILA, 
University of Colorado, for its hospitality during the completion of part of this work. 



- 19- 

REFERENCES 

Aharonian, F. 2000, astro-ph/0003159 

Antonucci, R., Hurt, T., & Kinney, A. 1994, Nature, 371, 313 
Blandford, R.D., & Levinson, A. 1995, ApJ, 441, 79 

Blazejowski, M., Sikora, M., Moderski, R., Madejski, G.M. 2000, ApJ, 545, 107 (B2000) 

Dermer, CD. 1995, ApJ, 446, L63 

Dermer, CD. & Schlickeiser, R. 1993, ApJ, 416, 458 

Emmering, R.T., Blandford, R.D., & Shlosman, I. 1992, ApJ, 385, 460 

Fossati, G., Maraschi, L., Celotti, A., Comastri, A., & Ghisellini, G. 1998, MNRAS, 299, 
433 

Ghisellini, G., Celotti, A., Fossati, G., Maraschi, L., & Comastri, A. 1998, MNRAS, 301, 
451 

Ghisellini, G., & Madau, R 1996, MNRAS, 280, 67 

Inoue, S., & Takahara, F. 1996, ApJ, 463, 555 

Komissarov, S.S., & Falle, S.A.E.G. 1997, MNRAS, 288, 833 

Konigl, A. 1981, ApJ, 243, 700 

Krichbaum, T.R, et al. 1998, A&A, 329, 873 

Kubo, H., et al. 1998, ApJ, 504, 693 

Lawson, A.J., McHardy, I.M., & Marscher, A.P. 1999, MNRAS, 306, 247 



-20- 

Lobanov, A.P. 1998, A&A, 330, 79 

Mannheim, K. & Biermann, P.L. 1992, A&A, 253, L21 

Maraschi, L. & Rovetti, F. 1994, ApJ, 436, 79 

Mattox, J.R., et al. 1997, ApJ, 476, 692 

Moderski, R., Sikora, M., & Bulik, T. 2000, ApJ, 529, 151 

Murray, N., & Chiang, J. 1997, ApJ, 474, 91 

Miicke, A., & Protheroe, R.J. 2000, astro-ph/0003282 

Nicastro, F. 2000, ApJ, 530, L65 

Padovani, P. & Urry, CM. 1992, ApJ, 387, 449 

Proga, D., Stone, J.M., & Kallman, T.R. 2000, astro-ph/0005315 

Rachen, J.P. 2000, astro-ph/0003282 

Rachen, J.P., & Meszaros, P. 1998, Phys. Rev. D, 58, 123005 
Sari, R., Piran, T., & Narayan, R. 1998, ApJ, 497, L17 

Sikora, M. 1997, AIP Conference Proceedings (4th Compton Symposium), 410, 494 
Sikora, M., Begelman, M.C., & Rees, M.J. 1994, ApJ, 421, 153 
Sikora, M., & Madejski, G. 2000, ApJ, 534, 109 

Sikora, M., Madejski, C, Moderski, R., & Poutanen, J. 1997, ApJ, 484, 108 
Sikora, M., Sol, H., Begelman, M.C. & Madejski, CM. 1996, MNRAS, 280, 781 
von Montigny, C, et al. 1995, ApJ, 440, 525 



- 21 - 

Wagner, S.J., et al. 1995, ApJ, 454, L97 
Wehrle, A.E. et al. 1998, ApJ, 497, 178 



This manuscript was prepared with the AAS macros v5.0. 



-22- 



FIGURE CAPTIONS 

Fig. 1. — Evolution of electron energy distribution. The upper panels are for radiative 
electron energy loss rates oc 7^ and independent of distance [7(r) = {dj/dt')rad — const]; the 
lower panels are for 7 oc the left panels are for zero adiabatic energy losses, the right 

panels are for SD-expansion adiabatic losses. The models are computed assuming p — 2.4, 
7^ = 0.1 jc, and Arinj/r^ — 1. Evolution is followed from t /to — 1 up to t /to — 3, and the 
time step between the presented curves is t/to — 0.1. 

Fig. 2. — Flare profiles: (a) for different values of the injection distance range, Avinj/ro; 
(b) for different adiabaticities. A; (c) for different electron injection time profiles, Q{t'); ; 
(d) for different jet opening-angles, 9j. Models (a), (c), and (d) are computed for 7 ^ 7c 
(radiative regime) and models (b) are computed for 7 <S 7c (adiabatic regime). Unless 
marked differently, the parameters used are: Avinj/ro = 1; Q{t') = const; A = 2/3; 9j — 
1/15. All models are computed for F = 15. 

Fig. 3. — Flares produced from double-triangle injection time profile. Panel (a) shows 
dependence on the injection distance range Arinj — 0.3, 1.0, 6.0 per triangle, while panel (b) 
is for different observers located at angle 9j — 1/60, 1/30 from the jet axis. 

Fig. 4. — Flare lags: (a) for different synchrotron frequencies; (b) for different ERC frequen- 
cies; (d) for different radiation processes, at frequencies as marked on (c). The model has 
been computed for the following parameters: bulk Lorentz factor, F = 15; initial injection 
distance, ro = 6.0 x 10^^ cm; injection distance range, Armj/ro — 1; constant injection 
rate Q{t') — const, with 7^ = 10, ^max — 6.0 x 10^ and p — 2.2, and normahzation 
given by = rUeC^ j d'y — 3.5 x 10^^ ergs s~^; adiabaticity A = 1; magnetic field 
intensity, B' = 0.4(ro/r) Gauss; and energy density of the ambient diffuse radiation field, 
ud — 3.5 X 10~^ ergs cm~^. 
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Fig. 5. — Modeling the Feb-96 event in 3C 279. Panel (a) shows model fit to the time averaged 
outburst spectrum. All observational data are simultaneous and taken from Wehrle et al. 
(1998). Panel (b): model 7-ray flare, with the flux integrated over energies > 100 MeV, 
vs. model X-ray flare, with the flux taken at 2 keV. Parameters of the model are: F = 20; 
ro = 6.0 X 10^''' cm; Ari„j/ro = 1; Q{t') = const, with 7^ = 27, ^rnax = 6.5 x 10^ and p = 2.4, 
and Lg = 1.6 x 10^^ ergs s~^; B' = 0.53(ro/r) Gauss; and energy densities of the ambient 
diffuse radiation fields, ubel — 4.9 x 10~''(ro/r)^ ergs cm~^ and um — 1.0 x 10~^ ergs cm~^. 
Panel (c) presents X-ray flare profile from panel (b) fitted to 2-lOkeV RXTE data from the 
campaign (Lawson, McHardy & Marscher 1999). 

Fig. 6. — Modeling the Jan-93 event in PKS 1406-076. Left panel: model spectra taken at 
the moment of the optical flux peak (dotted line) and at the moment of the 7-ray flux peak 
(solid line). The optical and 7-ray data are the fluxes averaged over the event. (Radio flux 
comes from a different epoch.) Right panel: model 7-ray flare, with the flux integrated over 
energies > 100 MeV. vs. model optical flare, taken in the R-band. The model parameters are: 
r = 15; ro = 1.2 x 10^^ cm; A^^./ro = 1; Q{t') oc (r - ro), with 7^ = 50, -^max = 2.2 x 10^ 
p = 2.2, and Lf,{t = 2io) = 3.5 x 10^^ ergs s"^; B' = 0.14(ro/r) Gauss; and energy density of 
the ambient diffuse radiation fields, um — 3.0 x 10~^(r — ro)/ro ergs cm~^. 

Fig. 7. — Variabihty around the ERG peak. Left panels are for the peak produced by 
electrons with energies 7 ~ 7c, while right panels are for the peak determined by electrons 
with 7 ~ 7^. The spectra are shown with a time step 5t/tQ = 0.14 and are followed up to 
t — 3to- The flares are shown at frequencies: (c) v — O.li^c (sohd hue) and v — lOOi^c (dashed 
hne), where Vc — (16/9)r^7^i>'7R; (d) v — Q.lvm (sohd line) and v — IQQvm (dashed line), 
where Um = (16/9)r^7^z/7R. {d'j/dt')ERc is assumed to be distance independent. 

Fig. 8. — The same as Fig. 7, but for {d'y/dt')ERC 1/r^. 
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